Amplitude death and synchronized states in nonlinear time-delay systems coupled 

through mean-field diffusion 



Tanmoy Banerje^] and Debabrata Biswa^j] 
Department of Physics, University of Burdwan, Burdwan 713 104, West Bengal, India. 
(Dated: May 2, 2013; Received :to be included by reviewer) 

We explore and experimentally demonstrate the phenomena of amplitude death and other syn- 
chronized states in coupled time-delay hyperchaotic systems interacting through mean-field diffusion. 
Studies in the parameter space reveal that, with increasing coupling strength, the coupled systems 
make a transition from the unsynchronized state to amplitude death state via in-phase (complete) 
synchronized states. We employ Krasovskii-Lyapunov theory to derive the sufficient parametric con- 
dition of occurrence of complete synchronization. Also, we carry out a linear stability analysis that 
predicts the region of amplitude death in the parameter space. The Lyapunov exponent spectrum, 
correlation functions, and eigenvalue spectrum are used to establish the occurrence of in-phase (com- 
plete) synchronization and amplitude death. To corroborate the analytical and numerical results, 
we demonstrate all the scenarios in an electronic circuit experiment; the experimental time-series, 
phase-plane plots, and generalized autocorrelation function computed from the experimental time 
series data confirm the occurrence of in-phase (complete) synchronization and amplitude death in 
the coupled oscillators. 
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Coupled dynamical systems show a plethora of 
complex collective behaviors like synchronization, 
phase-locking, amplitude death, etc. Amplitude 
death (AD) occurs in coupled oscillators when 
they interact in such a way as to suppress each 
others oscillations and collectively go to a stable 
fixed point. This intriguing and emergent phe- 
nomenon is common in nature and has practical 
importance in suppression of unwanted oscilla- 
tions in laser systems, electronic circuits, metal 
cutting tools, etc. Realizing the importance of 
this phenomenon, different techniques of induc- 
ing AD have been investigated. Several coupling 
schemes are proposed for this purpose each hav- 
ing their merits and demerits. In this context, 
mean-field diffusive coupling is one of the impor- 
tant coupling schemes because it removes the con- 
straint of having parameter mismatch or time- 
delay coupling to obtain AD. Interestingly, al- 
most all the previous studies on AD consider that 
the interacting systems do not posses an intrin- 
sic time-delay. Very few works are reported on 
AD in time-delayed systems, and most of them 
consider a delayed coupling scheme. Dynamical 
systems having intrinsic time-delay are infinite di- 
mensional and very complex, and thus they have 
to be treated separately. Further, practical imple- 
mentation of time-delay systems are difficult and 
challenging, thus studies on amplitude death in 
coupled time-delay systems and its demonstration 
in experimental electronic circuit are of consid- 
erable importance. In this paper, we extensively 
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study the dynamical behavior of time-delayed hy- 
perchaotic oscillators coupled with mean-field dif- 
fusion; we theoretically explore and experimen- 
tally demonstrate the phenomena of amplitude 
death and in-phase (complete) synchronization in 
the coupled systems. 



I. INTRODUCTION 

Cooperative behaviors in coupled dynamical systems 
are of significant interest in the field of physical science, 
biological science, and engineering applications [1 j. The 
prominent cooperative behaviors that occur in periodic 
and chaotic oscillating systems are synchronization pQ, 
amplitude death[2^, phase-flip bifurcation [3j, etc. Ampli- 
tude death (AD) is one of the fascinating and important 
emergent phenomena in which quenching of amplitude 
or cessation of oscillation to the steady state occurs in 
coupled dynamical systems under some proper paramet- 
ric conditions. Studies on AD has been attracting the 
attention of the researchers for more than two decades 
owing to its importance in the field of physical science, 
biology, oceanography, etc[2j[4]. The phenomenon of AD 
has first been reported by Yamaguchi and Shimizu [5]; 
later the same has been studied in detail by Bar-Eli [6], 
and Shiino and Frankowicz [7]. Recently, an extensive 
review on AD has been reported by Saxena et al. [2] that 
discussed several aspects of AD, and a thorough litera- 
ture review. These studies reveal that, depending upon 
type of the interacting systems and nature of the coupling 
schemes, two distinct bifurcation routes can give rise to 
AD: (i) Hopf bifurcation, where the origin of the un- 
coupled systems becomes stable beyond a threshold cou- 
pling strength, and (ii) saddle-node bifurcation, where 
new fixed points appear for sufficiently large coupling 
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strength. 

To induce AD in coupled oscillators, several coupling 
schemes have been proposed; a few of them are linear dif- 
fusive coupling [2 j, delayed coupling [8], dynamic coupling 
[2], environmental coupling [9 j, etc. In most of the cou- 
pling schemes, to induce AD, it is necessary that the sys- 
tems are either mismatched or having time-delayed con- 
nection. However, these constraints have been removed 
with dynamic coupling, conjugate coupling [10 , and lin- 
ear augmentation [TT] coupling, where AD occurs in iden- 
tical systems with instantaneous coupling. Recently, an- 
other interesting process of inducing AD in identical sys- 
tems has been reported by Sharma and Shrimali [12] that 
considers mean-field (MF) diffusive coupling. Earlier, 
amplitude death through MF coupling has been proposed 
in [13], and [14]; but those works were done for coupled 
limit cycle oscillators with either distributed frequencies 
or parameter mismatches. 

In all the above mentioned works on amplitude death, 
the coupled dynamical systems are considered to be low 
dimensional with no intrinsic time-delay. The presence 
of time-delay makes a dynamical system infinite dimen- 
sional and very complex [15]. Recent ongoing interest in 
time-delay systems originates from the fact that in the 
real world often we have to encounter with the systems 
having intrinsic time-delay; examples include, blood pro- 
duction in patients with leukemia (Mackey- Glass model) 
[16], dynamics of optical systems (e.g. Ikeda system) 
[TT] . population dynamics [18], El Niho/southern oscil- 
lation (ENSO) [19], etc. Due to infinite-dimensionality 
of the time-delay systems, studies on their collective be- 
haviors like synchronization and AD are much more in- 
volved, and need special attention [20] . In this context, 
synchronization of chaos and hyperchaos in time-delay 
systems is a well explored topic and extensive research 
has already been devoted. After the first report on syn- 
chronization of chaos in time-delay systems by Pyragas 
[21], several types of synchronizations in time-delay sys- 
tems have been reported; a few examples are: lag and 
anticipatory synchronization [22 , complete and gener- 
alized synchronization [23] , phase synchronization [24] , 
inverse-phase synchronization [25 , etc. 

Contrary to the phenomenon of synchronization, am- 
plitude death in time-delay systems is a less explored 
topic. The first observation of AD in time-delayed oscil- 
lators has been reported by Konishi et al. [26 in which 
dynamic and delayed couplings were studied. AD in net- 
works of delay-coupled delay oscillators is studied ana- 
lytically by Hofener et al. [27] . Refl28l reports AD in 
time-delayed oscillators coupled via multiple delay con- 
nections. Recently, an important technique of inducing 
AD has been proposed by Ghosh et al. [20], which em- 
ploy time-delayed open-plus-closed-loop coupling. In all 
of these works (except the case of dynamic coupling), 
AD is mediated by the presence of moderate or very long 
time-delay in the coupling path. The presence of time- 
delay in coupling path makes the dynamics of the cou- 
pled systems more complex, and at the same time much 



difficult for analysis and practical implementation. 

In this paper, our purpose is to study the amplitude 
death and other synchronized states in coupled identical 
nonlinear time-delay dynamical systems with mean-field 
diffusive coupling. The choice of the mean-field coupling 
scheme is motivated by the fact that, to induce AD in 
identical systems, it does not require a time-delay in the 
coupling path. In this category, as we have discussed ear- 
lier, the other popular schemes are conjugate coupling, 
linear augmentation and dynamic coupling; but conju- 
gate coupling is not possible in a first-order time delay 
system, and for other schemes, we require additional lin- 
ear system(s) to realize the proper coupling that, in turn, 
increases the dimension of the original system, and also, 
makes the coupling scheme complex for the experimental 
implementation. Further, the idea of mean-field coupling 
is widely used in physics (e.g., in three dimensional Ising 
spin model [29] interaction with mean-filed explains ferro- 
magnetism), laser applications [30 , biology [31], etc. 

In the present paper, we consider two identical non- 
linear first-order time-delay systems interacting through 
mean-field diffusion. We show that, depending upon the 
coupling strength and mean-field parameters, the coupled 
systems show a transition from the unsynchronized state 
to amplitude death state via in-phase and complete syn- 
chronized states. It is observed that the transition from 
the unsynchronized state to synchronized state solely de- 
pends upon the coupling strength and does not depend 
upon the mean-field parameter - mean-field parameter 
only controls the inherent dynamics of the coupled sys- 
tems. Transitions from the unsynchronized state to com- 
plete synchronized state is predicted analytically using 
Krasovskii-Lyapunov theory. Also, we carry out a lin- 
ear stability analysis on the synchronized state to find 
out the parameter zone where amplitude death occurs. 
To exemplify our study, we consider a prototype time- 
delay hyperchaotic system recently proposed in Ref. [32] 
The choice of this system is motivated by the fact that 
it shows hyperchaos even for a small or moderate time- 
delay, and it can be implemented with off-the-shelf elec- 
tronic circuit components. We compute the Lyapunov 
exponent (LE) spectrum of the coupled systems, correla- 
tion of probability of recurrence, cross correlation func- 
tions, and eigenvalue spectrum to identify the synchro- 
nized states and AD state in the parameter space. Using 
the eigenvalue spectrum of the coupled systems we show 
that the route to AD, in the present case, is through Hopf 
bifurcation. Finally, we set an electronic circuit experi- 
ment to demonstrate all the transition scenarios. Experi- 
mental waveforms and phase-plane plots are used to visu- 
alize the transition from the unsynchronized state to AD 
state via in-phase (complete) synchronized states; also, to 
confirm the occurrence of in-phase synchronization in the 
experimental circuit, we compute generalized autocorre- 
lation function using the experimental time-series data. 
We show that analytical and numerical results agree well 
with the experimental observations. 

The rest of this paper is organized in the following 
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way: The next section describes the mathematical model 
of the coupled time-delay systems with mean-field diffu- 
sion. Stability analysis is reported in Section |III| Nu- 



merical simulation results are described in Section IIVI 
Section |V| reports the experimental implementation of 
the coupled system and experimental demonstration of 
amplitude death. Finally, we summarize the main obser- 



vations of our study in Section VI 



II. MEAN-FIELD COUPLING 

We consider N number of first-order time-delay dy- 
namical systems interacting through mean-field diffusive 
coupling; mathematical model of the coupled system is 
given by 



= h(xi,x iTi ;p) + e QX - xi 



(1) 



with i = 1 • • • N, X = Y^iLi x i 1S the mean-field of the 
coupled system. x T = x(t — r), r £ R + is the constant 
time-delay, and p represents the m dimensional parame- 
ter space. The coupling strength is given by e, and Q is a 
control parameter that determines the density of mean- 
field p~2j [31] . Here the function h(xi,Xi Ti ;p) is given 
by h(xi,X2 Ti iP) — —axi — bif(xi Ti ), thus the individual 
units are represented by the following scalar first-order, 
retarded type delay differential equations: 



Xi -axi - bif(x iTi ), 



(2) 



where a > and bi are the system parameters, and 
is the intrinsic time delay associated with the individual 
systems. Eq.Q represents a general class of first-order, 
nonlinear, retarded delay-differential equations. For ex- 
ample, for Mackey- Glass system [16]: f(xi Ti ) = — 1+ ^ c ; 

for Ikeda system [T7] : f(xi Ti ) = sin(a^ r .), etc. Thus, 
Eq. represents mean- field diffusive coupling scheme for 
any first-order delay dynamical systems. 



III. STABILITY ANALYSIS 

In this section we analyze the asymptotic stability of 
the synchronization of the coupled systems given in Eq. 
([!]). Here we restrict our study to a pair (N = 2) of 
time-delay systems. 



A. Krasovskii— Lyapunov theory: complete 
synchronization 



The synchronization manifold of complete synchro- 
nization is x\ = x 2 . Let us define the error function as 
A = (x\ — x 2 \ and also let t\ — r 2 = r. Time evolution 
of the error function that describes the error dynamics of 
([I]) is given by 

A = -(a + e)A - (h - b 2 )f(x lr ) - b 2 f(x lT )A T . (3) 



This is an inhomogeneous equation and difficult to deal 
with; to make it homogeneous we impose the following 
constraint: 



bi = b 2 = b, 



(4) 



which is the necessary condition for complete synchro- 
nization. Now, Eq. ([3| becomes 



A = -(a + e)A-6/'(a;i T )A T . 



(5) 



According to the Krasovskii-Lyapunov theory [33] , a sta- 
ble synchronization implies the stability of the origin of 
([5]). The sufficient condition for the stability of synchro- 
nization requires the definition of a positive definite func- 
tional, V(i), given by 



V(t) 



A 2 + nj A 2 (t + ^)d^. (6) 



Here ji > is an arbitrary positive parameter. The sta- 
bility of the origin of (|J| requires that the time derivative 
of V(t) be negative. Now, 



dV 
~dt 



AA + /iA 2 - fiA 2 r1 

= -/iA 2 r(x,/i), 



A ' 



(7) 

where T(X, fi) = + b J^lA X + X 2 and X = ^ 

Thus, from Eq.Q it may be noted that the negativ- 
ity of dV/dt requires the following condition to be valid: 
Tmin > 0. Now, T min is derived as 

_ 4Ma + e-M)-fr 2 / /2 (*iJ 

J- rain — A o • \&) 



Hence T min > implies that 
b 2 f 2 (x lr ) 



e > 



4 M 



+ /i = $(/i). 



(9) 



Here 3>(/i) is a function of \i. Now, we find the minimum 
value $ min by setting || = 0. That gives \i = IM^Oi. 
With this value of /i, one gets the minimum value of $ 

as 



\bf'(xi r ) 



(10) 



Thus, Eq.(10) in Eq.rt9| gives the following sufficient con- 



dition of complete synchronization: 
a + e>\bf(x lT )\. 



(11) 



Note that, Eq.(ll) represents the sufficient condition of 



complete synchronization for any general first-order time- 
delay systems of the form given by Eq.Q coupled via 
mean- field diffusion. From Eq. (11) it is interesting to 



note that the mean-field parameter Q does not have any 
influence on the transition to the complete synchronized 
state, i.e., by varying Q alone, one cannot desynchronize 
the systems from the complete synchronized state. 
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B. Linear stability analysis: amplitude death 



IV. NUMERICAL SIMULATION 



Next, to find out the condition of amplitude death we 
analyze the stability of synchronization by considering 
the deviations from the synchronized state. The same 
for the low-dimensional systems (without intrinsic time 
delay) has been reported in Ref.9 and ReffT2l Let us de- 
fine ip and (j) to be the deviations from the synchronized 
states of the system variables x\ and x 2 in 0, respec- 
tively Then the linearization of the system along these 
deviations gives 

jJj = -aiP - &i/'(xi T )Vv + e(Q*±£ - (12a) 

=-a^-b 2 f\x 2T )^ T + e(Q^-cf). (12b) 

An exact analysis of Eq.([T2|) is not possible due to the 
presence of the delay term, which makes the charac- 
teristic equation a quasi-polynomial one. We consider 
b\ = 62, and for the the complete synchronization we have 
x\ = X2 and X\ T = x 2r = x T . Let us define g(x T ,i/j, ip T ) = 
-ai/j - bf(x T )i/j T , and g(x T , <j>, <j> T ) = -a<j> - bf r (x T )(j) T . 
Now, Eq.(|12|) reduces to 

», (13a) 
-0). (13b) 



ip = g{x T ,ijj,ijj T ) + e(Q 



^+0 
2 

= s(z r ,0,0 r ) + e(Q^ 



The Jacobian matrix of the system is described by 



f + e(g-l) 



6 + eQ-l) 



0, 



(14) 



where, we consider that the time- averaged values of 
g'{x T , ip, ^; r ) and g'(x T: 0, r ) are approximately same 
and are equal to an effective constant (5. This type of 
approximation has been used in Refs|T2j |34| for the low- 
dimensional systems without intrinsic time-delay; here 
we extend the same for the time-delay systems. 

Now the characteristic equation of the Jacobian matrix 



(14) is 



A - 2AA + A 



1 



= 0, (15) 
Thus, we have the following 

1), 



where, A = S + e( 
two eigenvalues: 

Ai =(5 + e(Q 
A 2 = S — e. 

For amplitude death to occur, Ai^ 
negative [12], which gives 



(16a) 
(16b) 
should be 



e 



(17a) 

e > d. (17b) 

Thus the critical parametric condition for which ampli- 
tude death occurs is given by 

Q cr = l , (18) 

along with e cr > 5; here Q cr and e cr are the critical 
values of the mean-field parameter and coupling strength, 
respectively. 



A. System description 

For numerical verification of the analytical predictions 
and demonstration of the collective behaviors, we con- 
sider the following first-order nonlinear retarded time- 
delayed system recently proposed in Refl32l 



x = —ax — bf(x T ), 



(19) 



where a and b are positive parameters. The nonlinear 
function f(x T ) is given by 

f(x T ) = —0.5n(\x T \ + x T ) + mtanh (lx T ), (20) 

where, n, m, and I are positive parameters that deter- 
mine the nature of the nonlinearity. There exists a large 
number of choices of these parameters for which chaos 
and hyperchaos can be observed [32] . 

The detailed chaotic and hyperchaotic dynamical be- 
haviors have been reported in Ref. 32. The system (19) 
(with ( |20| )) has only one (trivial) fixed point at x* = 0. 
It has been shown that, keeping b fixed, if one varies r, 
the system shows a period doubling route to chaos and 
hyperchaos (parameters are: a = 1, in = 1.15, m = 0.97, 
I = 2.19). For example, for b = 2.4, at r w 0.59, the 
fixed point loses its stability through Hopf bifurcation 
and a stable limit cycle appears; after a period dou- 
bling sequence, chaos and hyperchaos are observed at 
r w 1.61 and r w 2.54, respectively. The equivalent 
period-doubling route has been observed if one keeps r 
fixed and varies b; e.g., for r = 3, at b « 0.92 oscillation 
set in; chaos occurs for b ~ 1.71 and hyperchaos occurs 
for b w 2. Figure [TJa) shows the hyperchaotic attractor 
for b = 2.4 and r = 3. Figure [TJb) shows the Lyapunov 
exponent spectrum (along with the Kaplan- York dimen- 
sion) in the b parameter space for r = 3; the presence 
of multiple positive Lyapunov exponents along with the 
strange attractor ensures the occurrence of hyperchaos in 
the system. 



Numerical results 



The system equation ([!]) (with (19) and (|20|)) is sim- 
ulated numerically using Runge-Kutta algorithm with 
step size h = 0.01. The following initial functions 
have been used for all the numerical simulations: for 
the Xi-system: £ x i(t) — 0.95, and for the ^-system: 
€x2(t) = 0.85. Also, the following system design pa- 
rameters are chosen throughout the numerical simula- 
tions: a = 1, n = 1.19, m = 0.97, I = 2.19. We choose 
bi=b 2 =2.4, and r = 3, to ensure that both the systems 
be in the hyperchaotic region [32] . The numerical results 
are presented after discarding a large number of itera- 
tions to avoid transients. 

To observe the effect of coupling strength on the dy- 
namics of the coupled system, we vary e by keeping the 
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FIG. 2. Q=0.6: (a) Lyapunov exponent (LE) spectrum of 
largest five LEs. Inset shows the same, but now the range is 
e G (0.3,1.2)), and curves are smooth (through averaging) for 
clarity, (b) Bifurcation diagram of x\ with e. AD represents 
amplitude death region in e parameter space. 



FIG. 1. (a) Hyperchaotic attractor for b = 2.4 and r = 3. 
(b) Lyapunov exponent spectrum with b for r = 3, the vertical 
line is for b = 2.4; the inset shows Kaplan- York dimension 
(Dry)- Other parameter values are: o = 1, n = 1.15, m = 
0.97, I = 2.19. 



mean- field parameter Q fixed at Q = 0.6. At first, we 
compute the Lyapunov exponent (LE) spectrum of the 
coupled systems. Figure |2|a) shows the first five LEs 
in the e parameter space, e G (0,8). Inset of the figure 
shows the same in e G (0.3, 1.2), but here, for clarity, we 
use averaging of the data points to plot smooth curves. It 
can be seen that LE4 becomes negative at e ~ 0.65 that 
indicates the onset of phase synchronization pQ. Further, 
LE3 makes a transition from a positive value to negative 
value at e « 1, indicating the onset of complete synchro- 
nization. The transition of LE2 from a positive value to 
a negative value via null values indicates that increasing 
coupling strength makes the dimension of the coupled 
systems lower. With further increase in e, the largest 
Lyapunov exponent, LEI, becomes zero at e ~ 1.85, in- 
dicating the fact that the dynamics of the coupled sys- 
tems now become periodic. This periodic (and complete 
synchronized) state prevails until LEI becomes negative 
at e ~ 5.85. The transition of LEI from zero to a nega- 
tive value is indicative of amplitude death in the coupled 
systems. With further increase in e, LEI monotonically 




-0.2 -O.I 



(b) 



-0.3 



-0.2 -0.1 



FIG. 3. The eigenvalue spectrum of the coupled systems 
(a) e = 5 (i.e. near but before AD) (b) e = 7 (i.e., near but 
after AD); note that the real part of all the eigenvalues now 
become negative. 
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decreases toward a more negative value ensuring the sta- 
bility of the AD state. The AD state can best be observed 
from the bifurcation diagram of x x (Fig. [2jb)) with e, 
which is obtained by plotting the local maxima of x\. It 
is interesting to note, both from the LE spectrum and 
the bifurcation diagram, the presence of periodic win- 
dow (near e = 0.93), where LEI suddenly jumps to a 
zero value. 

Next, we compute the eigenvalue spectrum of the cou- 
pled systems to confirm the occurrence of AD. For this 
purpose we use the bifurcation package DDE-BIFTOOL 
[35] . Figures [3^ a) and (b) show the eigenvalue spectrum 
of the coupled systems at e = 5 (i.e., near but before 
AD) and e = 7 (i.e., near but after AD), respectively. 
It can be seen that, with increase in e, the real part of 
the largest eigenvalue changes from positive to negative 
value confirming the occurrence of AD in the coupled 
system. Further, we observe that, with increase in e, the 
largest complex conjugate pair of eigenvalues cross the 
imaginary axis 9(A) from right to left confirming that 
the route to AD is through Hopf bifurcation. 

The transition from the unsynchronized state to com- 
plete synchronized state through in-phase synchronized 
state is verified using correlation of probability of re- 
currence (CPR), and cross correlation function (CC). 
CPR is a quantitative measure of phase synchronization 
(PS) introduced in Ref. [36] and extensively used in the 
literature [24]. It is related with the generalized autocor- 
relation function (P(t)), which is defined as 



p( t ) 



, JVi-t 



WXi - x 



i+t\ 



(21) 



here 6 is the Heaviside function, Xi is the ith data point 
in the X variable, N\ is the total number of data points, 
e t is a preassigned threshold value, and ||.|| represents the 
Euclidean norm. CPR is defined as 36 



CPR: 



0~lO~2 



(22) 



where, P\^ 2 present that the mean value has been sub- 
tracted, and a 1)2 are the standard deviations of the Pi(t) 
and P2(t), respectively. For PS states, CPR w 1. Fur- 
ther, CC is defined as 



CC 



(( Xl (t) - ( Xl (t)))(x 2 (t) - (x 2 (t)))) 
^((x 1 (t)-(x 1 (t))f)((x 2 (t)-(x 2 (t))) 2 ) 



(23) 



CC is a measure of complete synchronization (CS)[T]; 
in the CS state, CC= 1. Figure [4] shows the variation 
of CPR and CC with e. Increase of both the measures 
from a zero value with increase in the coupling strength 
e agrees with the LE spectrum and bifurcation diagram. 
For e > 0.65, CPR attains values nearly equal to one 
indicating the onset of PS, and for e > 1, CC attains a 
value of unity indicating the onset of CS in the coupled 
systems. 
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FIG. 4. Variation of CPR and CC with e (Q=0.6). 



Finally, for the demonstration purpose, we show the 
time evolution of the system variables x\(t) and X2(t) 
(given in Fig {5] (middle column)), synchronization mani- 
fold in x\ — X2 plane (right column of Fig{5|, and phase 
plane plots in x±(t) — x\{t — r) plane (left column of 
Fig{5|, for different e values (with Q = 0.6). FigjEJa) 
shows that, for e = 0.2, the coupled systems are un- 
synchronized in the hyperchaotic zone. For e = 0.9 
(FigJEJb)), the systems are in-phase synchronized, and 
for e = 1.1 (FigjEJc)), the systems are completely syn- 
chronized. Further increase in e makes both the systems 
periodic (period-2 for e = 2 (d) and period-1 for e = 5 
(e)). AD is shown in Figj5^f) for e = 6.2, which shows 
that both of the coupled systems attain the zero steady 
state which is the only and trivial steady state of the 
uncoupled systems. 

Next, we keep e constant at e = 4, and vary Q, 
Q G (0, 1). We observe that, at this value of e, the cou- 
pled systems are always in the CS state. The variation 
of Q only controls the intrinsic dynamics of the coupled 
systems (in the CS state). To grab the system dynamics 
with the variation of Q, we plot bifurcation diagram of 
x\ with Q. It shows that, with increase in the mean- field 
parameter Q, the coupled systems make a transition from 
the AD state to chaotic and hyperchaotic state through 
a period doubling route. Figures {6^ a), (b), and (c) show 
the real time and phase plane plots for three representa- 
tive values of Q depicting the corresponding transitions. 
Figure j6^a) shows AD state for Q = 0.4, Figsj6^b) and 
(c) show the periodic and hyperchaotic oscillations for 
Q — 0.7 and Q = 0.9, respectively; also x\(t) —X2(t) plot 
in the right panel of these two figures indicate the CS 
state. 

In order to explore the dynamics of the system in the 
global parameter space we plot phase diagram in the 
e — Q parameter space. The LE spectrum, CPR, CC and 
similarity function [15] are simultaneously used to iden- 
tify different dynamical zones in the parameter space. 
Figure [7] depicts the phase diagram, which shows three 
distinct regions in the e — Q space, namely unsynchro- 
nized state (NS), in-phase or complete synchronized state 
(PS/CS), and amplitude death (AD) state. It is notewor- 
thy that the transition from the unsynchronized state 
(NS) to synchronized state (PS/CS) does not depend 
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x,(0-X](M) 



xi(t), x 2 (t) - 1 



Xi(t)-X 2 (t) 




FIG. 5. Q=0.6: variation of system dynamics and synchro- 
nization states for variable e. It shows that with increasing 
e the coupled systems make a transition from the unsynchro- 
nized state (hyperchaotic) (a) to amplitude death state (f) 
via hyperchaotic in-phase synchronized state (b) and hyper- 
chaotic or chaotic and periodic complete synchronized states 
(c.d.e). 



upon the mean- field parameter Q, but depends only upon 
the coupling strength, e, which is in accordance with the 
analytical result ( pTj ). Further, we plot the critical curve 
(dark solid (blue) line in Fig. [7]) for the transition from CS 
state to AD state using the analytical result ( 18 ) with an 
effective choice of S 



2.51; it matches exactly with the 
numerically obtained critical values in the phase diagram. 
The effective value of S is obtained by fitting numerical 
results with the analytical result (18). Next, in the phase 



diagram we show the threshold curve (dotted vertical 
line) for the transition from NS or PS to CS using analyt- 
ical result obtained in Eq.(ll); since f f (x\ T ) (of Eq.(ll)) 
is a time varying function, thus, for a general result, here 
we consider the maximum value of f'{xi T ). In the present 
case, \ f'(x\ T )\ = |— 0.5n(sgn(xi r ) + 1) + ml sech 2 (/xi x )|; 
in order to cope with the sgn() function, we consider 
the positive and negative values of xi T , separately. For 
x\ T > 0, we have \f'(x\ T )\ max = 1.14; for this value, 
we show the transition threshold line (dotted vertical 
line) in the phase diagram. It lies well within the nu- 
merically obtained region of CS. For x\ T < 0, we have 
\f (%i T )\max = 2.119; for this value we have a vertical 
line (not shown in the figure) at e = 4.08 that is also sit- 
uated well inside the complete synchronized zone in the 



1.6 
1 .4 

1 .2 

< 1 

0.8 
0.6 
0.4 
0.2 



8=4 



(b) Q=0.7 



1 .2 J 



(a) Q=0.4 



AD 




Q 



FIG. 6. e—4: Bifurcation diagram of x\ with Q. Insets show 
the time evolution of x\ and X2 (middle panel), x\ — x\ T (left 
panel), and x\ — X2 (right panel), (a) AD for Q = 0.4 (b) 
Periodic for Q = 0.7 (c) Hyperchaotic for Q = 0.9. 




FIG. 7. (Color online) Phase diagram in e—Q parameter 
space. AD: amplitude death; PS/CS: in-phase or complete 
synchronized state; NS: unsynchronized state. Solid line in- 
dicates the analytically obtained critical curve for obtaining 
AD (with 5 — 2.51), dotted vertical line indicates the syn- 
chronization transition curve. 



phase diagram. Since the Krasovskii-Lyapunov theory 
gives only the sufficient condition of stability of CS thus 
an exact prediction of the parameter values for the tran- 
sition to CS is not possible from Eq.(ll). Nevertheless, 



using this theory we can get a region in the parameter 
space where CS occurs. 



V. EXPERIMENT 

A. Electronic Circuit Implementation 

We set up an electronic circuit level experiment to im- 
plement the time-delay system (19) (with (20)) under the 
mean- field diffusive coupling scheme given by Eq. ([T]). 
Figure [8] shows the representative diagram of the exper- 
imental electronic circuit. The proposed circuit consists 
of three distinct parts, namely, the Xi-system (upper por- 
tion), X2-system (lower portion), and the circuit to real- 
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DELAY(T) 



FIG. 8. Representative diagram of the experimental circuit 
(see text for a detailed description). 



ize the mean-field coupling (middle portion). Both the x\ 
and ^-systems consist of a low-pass section (Rq — Co), 
nonlinear device (ND), delay block (DELAY), gain (bi 
and 62), and other circuitry used to realize the proper 
coupling. The ND block produces the nonlinearity of 
both x\ and ^-systems; the circuit to realize the ND 
block is shown in Fig.[9|a). For a given input voltage V T 
(say), this circuit has a nonlinearity, f(V T ), given by [32] 



f(Vr) 



\ + Vr) 



R5 
R 3 



PVsat tanh 



R\ Vsat 



(24) 



Here /3 and uo are certain scaling factors that depend upon 
the non ideal and asymmetric nature of the op-amps, and 
V sat is the saturation voltage of the op-amps. The gain 
part bi and 62 (= b) is realized with op-amp A3 as shown 
in the same figure. The delay part is implemented using a 
chain of cascaded active all-pass filters ( APF) [37 (shown 
in Figj9^b)) having a transfer function T(s) = — ai J=^, 
where, a\ is the flat gain of the delay block given by 
a\ = R12/R111 ^0 is the frequency at which the phase 
shift is 7r/2. Owing to the almost linear phase response, 
each delay block produces a time delay of To « RdC. 
Thus, n blocks of APFs produce an effective delay of 
Td = uRdC (n is an integer). 

The mean-field is realized by an inverting adder (A5) 
(of Figj8| with a gain equal to ~. The signals from 
Ro~Co nodes of x\- and a^-systems are fed to the invert- 
ing terminal of A5 via buffers. The gain of the inverting 




FIG. 9. (a) The nonlinear device (ND) along with b. (b) The 
delay block implemented by active all-pass filter (APF). (c) 
Experimental nonlinearity of the device, (for the parameter 
values see text). 



amplifier A6, (Rq/R), controls the mean- field parame- 
ter Q. Two difference amplifiers, A7 and A8, are used 
to realize the final stage of mean-field coupling. Note 
that, the coupling strength e is controlled by Riq/Rq. 
Buffers (BUF) are realized using unity gain noninverting 
amplifiers, and used whenever necessary for impedance 
matching and elimination of unwanted feedbacks. 

Let Vi(t) be the voltage drop across the capacitance 
Co of the low-pass section of xi-system, and that of x 2 - 
system be V 2 (t). Then the equations that represent the 
circuit dynamics are 



RqCq 



RqCq 



dVijt) 
dt 



dV 2 (t) 
dt 



Rio ( Rq 

Rg \ R 

Rio 

Rg \i?6 

^10 1 Rq 



-v 2 (t) 



V{t)-V l {t) 

R7 f(v> TD . 



(25a) 



Rq \ R 



V(t)-V 2 (t) .(25b) 



here f(V lTo ) = f(V z (t-T D )) (i = 1,2) is given by Eq. 



■ _ Vj(*)+Vb(*) 



24 



and V(t) - 

Now we define the following dimensionless parameters 



and variables: t - 

Vl r„ y 2 (t) 



RnCn 



_ Rq_ 



v sc 

h = u>% , b 



V sa - 
Re' 1 



R0C0' Xl _^V sat ' 

x - ^ 



Rx_ 

R 8 



Vsat 



m 

1, e 



R* 

R 4 ' 

Rio 

Rg ' 



and Q = With these, (25) reduces to the following 
dimensionless form: 



dx\ 
~dtt 
dx2 
~dt 



- Xl (i) - bf(x u ) + e(QX - Xl ), (26a) 
-x 2 (t) - bf(x2 T ) + e(QX - x 2 ), (26b) 



with 



f(Vr) 



-0.5ni(\v T I 



- mi tanh (l\v T ), (27) 
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where, v = xi : X2- Thus, Eq.(|26|) (with (27)) is equivalent 
to Eq.Q (with [19] and ([2Q|)) with a = 1 and proper choice 
of Tii, and Zi. 



B. Experimental results 

Since we are considering the identical x\- and #2- 
systems, thus the resistance values in both the circuits 
are kept same: R x = 10 kft, i? 2 = 18.55 kft, i? 3 = 18.55 
kft, i? 4 5.6 kft, R 5 10 kft, i? 6 = 1 kft, i? 8 = i?i = 1 
kft. In the coupling part of FigJSJ i? = 10 kft. The low- 
pass sections have i?o = 1 and Co = 0.1 /iF. The APF 
section of Figgb) has R n = R 12 = 2.2 kft, C = 10 
nF, and i? D = 10 kft. All the op-amps are TL 074 IC 
(quad JFET op-amp) with ±15 volt power supply. The 
resistors and capacitors have 5% tolerance. Rg and Rq 
are varied with precession potentiometers (POT). 

With these values, the experimental nonlinearity is 
shown in Figj9^c), which is same for both the systems. 
To drive the systems into hyperchaotic zone we set r = 3 
(using three stages of APFs), and 61=62=2.4 by setting 
R 7 = 2.4 kft. 

(i) First we consider the case of variable e and fixed 
Q. We set Rq = 7.5 kft and vary e by varying Rg. 
The results of this variation are shown in Fig. [TUJ where 
the left column shows the phase plane of the xi-system 
(in V\{t) — Vi(t — Td) plane), the middle column shows 
the waveforms of V\(t) and V2(i), and the right column 
shows the phase plane plot in V\(t) — V^(t) plane, (a) For 
Rg = 10.77 kft, the scenario is shown in the first row a(l- 
3) of FigjlOj (al) shows the hyperchaotic attractor, and 
(a2) and (a3) show that there is no correlation between 
the coupled systems for these parameter values, and both 
the systems evolve independently, (b) For Rg = 862 ft, 
the behavior is shown in the second row, b(l-3); (b2) 
and (b3) suggest the occurrence of in-phase synchroniza- 
tion, (c) The third row (c(l-3)) shows the situation for 
Rg = 120 ft; (c2) and (c3) show that both the systems 
are in complete synchrony with each other. (d,e) Fur- 
ther decrease in Rg, i.e., increase in e, causes the systems 
to enter into the state of periodic oscillation. Period-2 
(fourth row d(l-3)) and Period- 1 oscillations (fifth row 
e(l-3)) are shown for Rg = 86 ft and Rg = 30 ft, re- 
spectively. It can be observed that in periodic states also 
the systems are in complete synchrony, (f) At very low 
coupling resistance the coupled systems show amplitude 
death (AD); Fig 10^f2) shows the waveforms of the cou- 
pled systems for Rg = 15 Vt that indicates the occurrence 
of AD, i.e., now the oscillations in both the systems die 
out; for better visibility of both the waveforms we add an 
offset to V2(t) (lower trace) with respect to V\(t). The 
phase plane plots FigjlO^fl) and (f3) show the collapse 
of oscillation to the steady fixed point, confirming AD. 
Thus, we have demonstrated that, fixing Q to a high 
value (< 1), when coupling strength is varied from low 
to high values, the systems show a transition from an 
unsynchronized state to the amplitude death state via 





(e2) 



<f2> 




FIG. 10. (Color online) Experimental waveforms and phase- 
plane plots with fixed Q (Rq = 7.5 kQ) and variable e. The 
left column shows the phase-plane plots (Vi(t) —Vi(t — Td)), 
the middle column shows the real time waveforms of Vi(t)- 
(yellow) and V2(t)-(blue), and the right column shows the 
phase-plane plots in (Vi(t) — V^(t)) plane. a(l-3) The un- 
synchronized state for Rq = 10.77 kf2; b(l-3) in-phase syn- 
chronizatin for Rg — 862 Q; c(l-3) complete synchronization 
for Rg = 120 ft; d(l-3) period-2 oscillation at Rg = 86 ft; 
e(l-3) period- 1 oscillation at Rg = 30 ft, and f(l-3) show the 
amplitude death for Rg = 15 ft; in (£2) the trace of 
(blue) is shifted downwards by 1.2 volt from that of V\{t) 
(yellow trace). (Scale div: Second column (a2-f2): x-axis: 25 
/isec/div, y-axis:1.25 volt/div. Other plots: x and y-axes: 0.5 
v/div) 



in-phase and complete synchronized states. 

For the confirmation of occurrence of in-phase synchro- 
nization we compute P\(t) and defined by Eq.(|2T|) 
from the experimental time-series data (acquired using 
DSO, Tektronix TDS2002B, 60 MHz, 1 GS/s). We 
choose the threshold value e t = 0.01, and Ni = 2400. 
Figure pJ^a) shows P(t)s for the unsynchronized case 
(with the same parameter values used in Fig. 10 'al-a3)); 
from the figure we can see that peaks of Pi (T) does not 
match with that of in the t-axis, confirming the 
fact that the phases of the two oscillators are not syn- 
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FIG. 11. (Color online) Plots of generalized autocorrelation 
functions using the experimental time-series data (a) unsyn- 
chronized state (parameters are the same as Fig ; 
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al-a3)), 

(b) in-phase synchronized state (parameters are the same as 
Fig[Mbl-b3)). 



chronized. Figure [TTJb) is for the parameter values used 
in Fig. 10 'bl-b3)); here we see that the dominant peaks of 



P\{t) and P2(t) matches exactly in the t-axis, confirming 
the in-phase synchronization of the coupled oscillators. 

(ii) Next, we consider the case of fixed high e and a 
variable Q. For this purpose, we set R$ = 139 and 
vary Rq from a low value to high values. The whole sce- 
nario is shown in F igfl2| (interpretation of three columns 
is the same as Fig|10|). (a) For a low value of Q, AD 
occurs; the situation is shown in the first row (i.e., a(l- 
3)) of Fig. " 
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for Rq = 8.42 kQ. It can be seen that 
both the systems are in the AD state. In (a2), an offset 
is kept between the trace of V\(t) and for clar- 

ity, (b) Increase in Q results in oscillation; this is shown 
in the second row of Fig. [l2] for Rq = 8.76 kft; (b2) 
and (b3) show that both the systems are in complete 
synchrony in period- 1 state, (c) Further increase in Q 
results in period doubling; the third row (c(l-3)) shows 
this at R Q = 9.52 kft. (d) With a large Q, both of the 
systems enter into a chaotic or hyperchaotic zone; Figs. 
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d(l-3)) show this for Rq = 9.76 k^. It is noteworthy 
that in all these cases both of the systems are in the com- 
plete synchronized state (as shown by the right column 
of Fig. 12) irrespective of the value of Q. Thus, we have 



demonstrated that for a sufficiently large value of e the 
variation of Q from a low value to a high value brings the 
system from the AD state to hyperchaotic state. 



VI. CONCLUSION 

We have explored the phenomenon of amplitude 
death in time-delay nonlinear dynamical systems cou- 
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FIG. 12. (Color online) Experimental plots for fixed e 
(R 9 = 139 Q) and variable Q. a(l-3) AD for Rq = 8.42 
kQ (the plot of V2 (blue) is shifted downwards by 1.25 volt 
from that of the Vi for clarity); b(l-3) period- 1 oscillation for 
R Q = 8.76 kQ; c(l-3) period-2 oscillation for R Q = 9.52 kQ; 
d(l-3) hyperchaotic oscillation for Rq =9.76 kQ. In all the 
cases the right column shows that the systems are in com plet e 
synchronized state. (Scale divisions are the same as Fig[lO| 



pled through mean-field diffusion. It has been shown 
that, depending upon the values of coupling strength and 
mean-field parameters, the coupled systems show a tran- 
sition to the amplitude death state from an unsynchro- 
nized state via in-phase (complete) synchronized states. 
A general stability analysis has been carried out using 
Krasovskii-Lyapunov theory to derive the sufficient con- 
dition of complete synchronization. Also, to predict the 
zone of amplitude death in the parameter space, we have 
performed a linear stability analysis around the synchro- 
nized states of the coupled systems. We have exemplified 
our results numerically using a prototype hyperchaotic 
time-delay system. Through the LE spectrum, correla- 
tion functions, and eigenvalue spectrum we have identi- 
fied the zone of amplitude death in the parameter space. 
It has been found that the numerical results agree well 
with the analytical derivations. The eigenvalue spectrum 
of the coupled systems revealed that, in the present sys- 
tem, the route to amplitude death is through Hopf bifur- 
cation. Finally, we set an experiment using electronic cir- 
cuit to demonstrate all the phenomena including ampli- 
tude death. It has been observed that the experimental 
results qualitatively agree well with the analytical results 
and numerical observations. 

In the present study, we have considered the identical 
systems; parameter mismatch in terms of the intrinsic 
time-delays or system parameters gives rise to a more 
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complex behavior in the coupled systems. Further, pos- 
sibility of amplitude death in an arbitrary number of 
coupled time-delay systems in a network under differ- 
ent types of coupling schemes needs further study, and 
studies are going on in that direction. 
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